EPJ manuscript No. 

(will be inserted by the editor) 



Coherence length of an elongated condensate: a study by 
matter-wave interferometry 



M. Hugbart la , J. A. Retter 1 , F. Gerbier 2 , A. Varon 1 , S. Richard 3 , J. H. Thywissen 4 , D. Clement 1 , R Bouyer 1 , A. Aspect 1 

1 Laboratoire Charles Fabry de lTnstitut d'Optique, UMR8501 du CNRS, 91403 Orsay Cedex, France 

2 Johannes Gutenberg-Universitiit, 55099 Mainz, Germany 

s THALES Research & Technology, 91404 Orsay Cedex, France 
4 Department of Physics, University of Toronto, Canada 

Received: date / Revised version: date 



Abstract. We measure the spatial correlation function of Bose-Einstein condensates in the cross-over region between 
phase-coherent and strongly phase-fluctuating condensates. We observe the continuous path from a gaussian-like shape 
to an exponential-like shape characteristic of one-dimensional phase-fluctuations. The width of the spatial correlation 
function as a function of the temperature shows that the condensate coherence length undergoes no sharp transition 
between these two regimes. 
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1 Introduction 

The high atomic phase-space density provided by a Bose- 
Einstein condensate (BEC) 1 1 1 has driven interest in guiding 
atoms in a manner analogous to guiding laser light through 
single-mode optical fibres. The most advanced technology to 
date is based on atom chips, where the fields which trap and 



is observed |27|: below a characteristic temperature deter- 
mined by the atom number and the trapping frequencies, the 
condensate is nearly phase coherent, but above the popu- 
lation of the axial modes is high and phase fluctuations may 
be pronounced. For weakly elongated condensates, can be 
higher than the transition temperature T c , so that the conden- 
sate is nearly phase coherent at all temperatures. In contrast, 
guide the atoms are created by microfabricated structures I2I3I4I5 6 Bbi chips can easily produce traps with high aspect ratios 



Many groups around the world have already succeeded in prepar- 
ing BEC on an atom chip 18 1911011111211311 41 and one can en- 
visage using this technology to create integrated atom inter- 
ferometers I15I16I17I18I . In this context, a precise characteri- 
zation of the phase coherence properties of the condensate is 
crucial II19I20I21I . 

Trapping quantum gases on atom chips naturally involves 
highly elongated, guide-like traps. Changing the dimensional- 
ity of the system from three-dimensional (3D) towards one- 
dimensional (ID) has a profound effect on the phase coherence 
of the condensate. In 3D condensates of modest aspect ratio, 
experimental results show that phase coherence extends across 
the whole cloud [22 23 1, even at finite temperature |24|. How- 
ever, in the ID regime, thermal excitations of the low energy 
axial modes lead to phase fluctuations which degrade the phase 
coherence |25 26 1. Condensates exhibiting such phase fluctua- 
tions are known as quasi-condensates. 

In the intermediate regime of elongated 3D condensates 
with a high aspect ratio, a behaviour similar to the ID case 
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(~ 1000) for which can be much smaller than T c . Phase 
fluctuations are therefore likely to impose limits on the perfor- 
mance of atom chip devices and need to be well understood. 

In an elongated condensate, the wavelength of the low en- 
ergy axial modes is longer than the radial size of the conden- 
sate, so these excitations have a ID character. However, the 
wavelength of these excitations can be much shorter than the 
axial length of the condensate, reducing the phase coherence 
length in the axial direction. An important feature of quasi- 
condensates is that density fluctuations remain suppressed in 
the trap, due to the mean field energy, even in the presence 
of large phase fluctuations 1 :271211 . Therefore, the 3D quasi- 
condensate has the usual parabolic profile in the Thomas-Fermi 
limit. 

Phase-fluctuating condensates in elongated traps were first 
observed by the conversion of phase fluctuations into density 
fluctuations after a sufficiently long free expansion 1 19 1 and by 
a condensate-focussing technique |20|. Quantitative measure- 
ments of the phase coherence length have since been obtained 
from the momentum distribution 1 2 1 1 and the second-order cor- 
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relation function 1281 . The results of each of these experiments 
showed good agreement with theory 1 27 1 in the strongly phase- 
fluctuating regime (at temperatures T ^> T^). However, neither 
experiment explored the cross-over region (T ~ X^) between 
phase-coherent and strongly phase-fluctuating condensates. 

In this article we describe a new experiment using a matter- 
wave interferometer and Fourier-space analysis to measure the 
spatial correlation function, thereby extending our measure- 
ments into the cross-over region (T ~ T^). Our results agree 
with the predicted shapes of the correlation function: for T ^> 
T^, we find exponential-like correlation functions as predicted 
for significant phase fluctuations, whereas at T ~ T$ we find 
a gaussian-like shape, as expected when the phase profile is 
nearly flat and the correlation function decay is dominated by 
the density profile. The coherence length as a function of T/T^ 
follows the trend predicted by theory, showing that the coher- 
ence length increases smoothly as the temperature falls and 
that there is no sharp transition at Tq. This highlights the fact 
that phase fluctuations occur at all finite temperatures, even if 
these effects are too small to be resolved experimentally for 
more spherical traps. However, whereas our previous measure- 
ments based on momentum spectroscopy 1211 . realized for high 
T/T,!,, were in full agrement with the theory, two observations 
remain unexplained in the interferometric method. First, as in 
a previous experiment 1 22 1, our experimental measurements of 
the coherence length are shifted from the theoretical prediction, 
by about 20% for T /T^ = 0, even after taking the limitations 
of our imaging system into account. Second, our interferome- 
ter produces unexplained supplementary fringes outside of the 
region where the condensates overlap, and we note that similar 
unexplained fringes appear in other published data 1 29 1 . These 
supplementary fringes do not seem to be compatible with in- 
terference of the thermal cloud observed in 1 30 1. 

2 Measurement of the coherence length by 
atom interferometry 

A natural method to study the coherence length along the long 
axis 2 of a condensate, or a quasi-condensate, is to use atom 
interferometry. With atomic beam-splitters, one produces two 
daughter copies of the initial condensate with a separation s, 
and observes the interference pattern appearing in the atomic 
density: 

n(r) cx |iPb(r - ^e z ) + e*+**9 (r + |e,)| 2 

(x|^ (r-^e z )| 2 + |^(r + ^e 2 )| 2 

+ 2Re[e i ^ *(r - S -e z )M>»{r + |e,)], (1) 

where e z is the axial unit vector, ^n(r) is the wavefunction de- 
scribing the initial condensate and <f> re \ a relative phase shift 
produced by the interferometer and the free fall of the conden- 
sate. 

Let us first consider the behaviour of a fully phase coherent 
condensate. During free expansion, it acquires a phase distri- 
bution proportional to z 2 1311 . The phase difference <fi le \ be- 
tween two displaced copies of the condensate is therefore pro- 
portional to zs, giving rise to an interference pattern of straight 



fringes, uniformly spaced along the longitudinal direction, the 
spatial frequency of the fringes being proportional to the sepa- 
ration s. The fringe contrast integrated over the entire conden- 
sate gives the first-order correlation function at s: 

CW(s) = /> r <^(r-|e z )<? (r+|e z ). (2) 

Repeating this contrast measurement for different separations 
s, one can study the decay of C^(s) with increasing s |22|. 
For a fully phase coherent condensate, the first-order correla- 
tion function decay reflects only the width of the density profile 
n(r) E|. 

In the case of a quasi-condensate, ID thermal excitations 
cause the phase to fluctuate along the longitudinal axis, both 
spatially and temporally. In our experiment, these fluctuations 
are small compared to the parabolic phase developed during 
free expansion. Therefore when we image the overlapping con- 
densates after free expansion, we still observe straight fringes, 
but they are no longer strictly periodic. Small local phase shifts 
add a "jitter" to the fringe spacing, which in Fourier space has 
the effect of broadening the peak at the spatial frequency of 
the fringes and thereby reducing its height. As s is increased, 
the fringes are perturbed more strongly, because the condensate 
phase becomes less correlated at larger separations. The con- 
trast therefore decreases faster with s than in the fully phase- 
coherent case. The greater the amplitude of the phase fluctu- 
ations, the faster the contrast decreases with s. Therefore by 
measuring the width of the correlation function at different 
temperatures, we extract the temperature dependence of the co- 
herence length. Further information is obtained from the shape 
of the correlation function 113 31 . 

In the presence of phase fluctuations, each realization of 
the experiment gives a different interference pattern, even with 
fixed experimental conditions. Expression (0 must therefore 
be generalized to: 

C«( S ) = J d 3 r (S£(r - ^e z )%(r + |e z )>, (3) 

where the brackets () denote a statistical average of the random 
process describing the random phase. In practice, one must re- 
peat the experiment at a given separation s and average the 
contrast measurements over many quasi-condensates to obtain 

This principle was used by Michelson in his famous astro- 
nomical interferometer, whose goal was to measure the spa- 
tial coherence of the light field arriving from a star, in order 
to deduce the diameter of the star |34|. However, Michelson's 
method is plagued by the existence of a randomly fluctuating 
relative phase between the two inputs of the interferometer, and 
various methods insensitive to the relative phase fluctuations 
had to be developed 1 35 1. A similar problem appears when one 
tries to determine the coherence length of a condensate or a 
quasi-condensate with interferometry. In addition to the con- 
trolled relative phase between the two components interfering 
in Q, there is an uncontrollable relative phase due to exper- 
imental problems such as a residual phase shift between the 
lasers creating the two copies or a random velocity kick im- 
parted to the sample. In order to overcome this problem |36|, 
D. Hellweg et al. |28| have used an analysis analogous to the 
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Hanbury Brown and Twiss method 1 37 1, since it is based on the 
measurement of the second-order correlation function which is 
insensitive to global phase shifts. In contrast, our method is 
in line with the initial method of Michelson who could visu- 
ally evaluate the contrast of the randomly moving fringes he 
was observing. The decrease of this contrast as a function of 
the telescopes' separation gave a direct measurement of the co- 
herence length. Similarly, we directly evaluate the contrast of 
the fringes by taking the modulus of the Fourier transform of 
the fringe pattern. The decrease of that contrast as a function 
of the separation s yields the coherence length of the quasi- 
condensate. 



3 Experiment 

3.1 Creation of elongated Bose-Einstein 
condensates 

In our experimental setup 1 38 1, a Zeeman-slowed atomic beam 
of 87 Rb is trapped in a MOT, and after optical pumping into 
the §S\/2\F = l,mp = —1) state is transferred to a mag- 
netic Ioffe-Pritchard trap created by an iron-core electromag- 
net. Our design allows us to lower the magnetic field at the 
bottom of the trap to a few Gauss and thus to obtain very tight 
radial confinement |39|. Using this trap, we are able to cre- 
ate condensates very close to the ID Thomas-Fermi regime 
1 40 1, as was demonstrated in [41 42 43 1 . In the present exper- 
iment, we produce condensates further into the 3D regime so 
that we can explore the cross-over regime (T ~ T$). We use 
two different trap configurations: in the first, the final radial 
and axial trap frequencies are respectively uo±_ — 2tt x 395 Hz 
and oj z = 2tt x 8.67 Hz, giving an aspect ratio of 45; for 
the second trap configuration, the final frequencies are ui± = 
2ir x 655 Hz and oj z = 2ir x 6.55 Hz, with aspect ratio 100. 
In this way, we obtain needle-shaped condensates containing 
around 3 x 10 5 atoms, with a typical half-length L ~ 85 /.tm 
in the first trap and L ~ 120 /im in the second. We control the 
final number of atoms by holding the condensate for a variable 
time, typically a few seconds, in the presence of an rf shied. The 
absolute number of atoms is calibrated from a measurement of 
the critical temperature, taking into account the effects of inter- 
actions 1441 . For condensates with small condensate fractions 
(less than 60 %), the temperature is obtained by fitting a gaus- 
sian distribution to the thermal wings of the cloud. The tem- 
perature is then extrapolated from the final frequency of the rf 
ramp to lower temperatures for which the thermal fraction is 
indiscernible I45I43I . 

3.2 Interferometry set-up and timing 

As shown in Figure ^ we implement the interferometer using 
a sequence of two -| -Bragg pulses, which act as matter-wave 
beamsplitters. The set-up consists of two laser beams counter- 
propagating along the longitudinal trap axis, each of intensity 
~ 2 mW cm" 2 , red-detuned by A = 6.6 GHz from the 87 Rb D2 
line at A = 780.02 nm. Two acousto-optic modulators driven 
by frequency synthesizers produce a small relative detuning 
S, tuned to the two-photon Bragg resonance S = 2hk1/m, 
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Fig. 1. The condensate is released from the trap at t = and has a half- 
length of L. A sequence of two Bragg pulses is applied to generate the 
interferometer. Two output ports are created (p — 0, p = with 
complementary fringe patterns. 



with rn the atomic mass and = 2tt/X. The momentum 
width along z of the expanding condensate corresponds to a 
frequency width of 200 Hz. Therefore we use Bragg pulses of 
~ 100 /is, short enough such that the corresponding 1.6 kHz 
frequency width (full width at half maximum) is sufficient to 
couple the entire condensate. The thermal cloud surrounding 
the condensate has a momentum distribution with a frequency 
width ranging from 12 kHz to 60 kHz, much larger than that of 
the condensate. Thus, only a small fraction of the thermal cloud 
is coupled by the Bragg pulses |45 1. By controlling the Bragg 
pulse length, we realize a ^ -pulse which splits the condensate 
into a coherent superposition of two wavepackets with veloci- 
ties differing by 2wr = 2hk\ J jm = 11.72 mmr 1 , where i>r 
is the 2-photon recoil velocity. The interferometer sequence is 
illustrated in Figure^ The condensate is held in the trap for at 
least 2 s at the end of the final rf evaporative-cooling ramp, to 
allow residual oscillations to be damped 12 11201 . After switch- 
ing off the trap, the condensate is allowed to expand freely for 
2 ms before the first 5 -pulse is applied. During this expansion 
the condensate density reduces by two orders of magnitude, 
so collisions between the diffracted wavepacket and the origi- 
nal condensate become negligible. During a free-evolution time 
2ms< T s < 10 ms, the two wavepackets separate to a distance 
s = 2v-rT s . The second \ -pulse completes the interferometer, 
and we observe interference in each of the two output ports, 
which differ in momentum by p = 2Kk\ J . The condensate is 
imaged by absorption perpendicular to the long axis z after a 
29ms total time-of-flight l46l . 

For a given set of experimental conditions (condensate atom 
number and temperature), the experimental correlation func- 
tion is acquired by taking a sequence of interference images 
with different condensate separations s ranging from 0.2L to 
1.2L, varied by changing T s . For smaller separations, we do 
not observe enough fringes to obtain a reliable measurement 
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of the contrast. At the maximum value of s, the contrast has 
reduced such that the fringes are no longer discernible above 
the noise. Typical images for 0.3L < s < 1.1L are shown 
in Figure 13 For each value of s, typically 5 images are taken, 
so that a statistical average can be performed. The fringe con- 
trast is then measured, giving the correlation function. Corre- 
lation functions have been obtained at various temperatures T 
between 100 and 230 nK and for condensate atom numbers No 
between 0.5 x 10 5 and 2.5 x 10 5 . These conditions correspond 
to 0.8 < T/T < 8, where T = 15h 2 N /16mk B TL 2 ED- 



4 Analysis of interferograms 
4.1 Interferogram 

As shown in 125 1 , a quasi-condensate is well-described by a 
fluctuating complex field ip(p,z) — \J n(p, z)e l ®( p,z \ with 
fixed density distribution n(p, z) and fluctuating phase <P(p, z). 
In the following, ^(p, z) represents the wavefunction of the 
condensate after the free-fall expansion. The first Bragg pulse 
is applied after 2 ms of free expansion, at which time the den- 
sity has reduced such that interactions between the atoms are 
negligible. We therefore assume that the different copies of 
the condensate propagate independently. The phase distribu- 
tion can be expressed as <P(p, z) = az 2 + f3p 2 + 4>th{z), where 
4>th{z) represents the thermal phase fluctuations and the quadra- 
tic terms represent the parabolic phase developed during expan- 
sion. 

We now consider the interference pattern produced at one 
of the output ports of the interferometer. For a separation s, we 
obtain the atomic density distribution: 

1 



s (urn) 



n out (p,z) = ^\V( Pl z - s/2) + V{p, z + s/2)\ 2 



71- 

~~4~ 



cos[A<P(z) 



(4) 



with n± = n(p, z ± s/2) and A$(z) = <P{z + s/2) - <P(z - 
s/2). The global phase shift <fi g is due to random, uncontrolled 
phase shifts between the two Bragg pulses. The phase differ- 
ence between the two copies is A<P(z) = azs + (f>th{z + s/2) — 
<pth(z — s /2)- The density of the condensate is small when the 
first Bragg pulse is applied, so we neglect a small relative ve- 
locity due to repulsion between the two copies. At the second 
output port, the two condensate copies have an additional it rel- 
ative phase shift due to the Bragg pulses, thereby producing a 
complementary fringe pattern. In our data analysis, the images 
of each output port are treated separately. 

Since the global phase shift <fi g fluctuates from shot to shot, 
we cannot average over different images at the same separation 
s. Instead, we take the contrast of each image individually and 
then average the contrast. 



4.2 Analysis in Fourier Space 

The atoms are imaged by absorption along the vertical y-axis, 
perpendicular to the long axis of the trap. The image we obtain, 
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Fig. 2. Left: One of the two output ports of the interferometer for 
different separations s between the two condensate copies. The half- 
length of the initial condensate is L = 85p m. Right: Profiles of the 
2D Fourier transform (absolute value) of each of the images on the 
left. The profiles are normalized to equal central peak height (propor- 
tional to total atom number). The position of the second peak (delim- 
ited by the two lines) gives the spatial frequency of the fringes, and 
the ratio of the height of the second peak to the central peak gives the 
fringe contrast. 



rescaled to units of 2D atomic density, is the integrated density: 
I(x,z) = const / dyn out . (5) 



We take the 2D Fourier transform of this image and extract its 
profile along the zero radial frequency k x — axis: 



m k z 



const / d rn ou te 



ik z z 



(6) 



Typical images and their 2D Fourier transform profiles are shown 
in Figure|2]for different separations s. The contrast of the fringe 
pattern is given by the ratio 2J[0, fco(s)]/J[0, 0], where fcrj(s) — 
as is the dominant spatial frequency of the fringe pattern. The 
profiles of Figure |2] show clearly the increasing spatial fre- 
quency and decreasing contrast as a function of s. 

To extract the correlation function, we take the complex 
amplitude of the Fourier peak at the spatial frequency fco of the 
fringes: 



l[0,k (s)} =e 1 ^ / d 3 r^npZe 



(7) 



In the absence of the global phase shifts <f> s , the correlation 
function (s) would be obtained by taking the statistical av- 
erage of equation 0: 



(I[0,ko(s)}) 



d 3 r {y/r. 



JA<P th \ 



(8) 



which is identical to relation l|3}- However, </> g changes from 
shot to shot and prevents us from averaging Fourier transforms 
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Fig. 3. Points represent the contrast extracted from our simulated ab- 
sorption images (see text) as a function of s/L, L being the conden- 
sate half-length, for T/T<t, = 1. These points are fitted by the product 
of a gaussian and an exponential (solid line) which we take to be our 
theoretical effective correlation function C cff (s). 



directly in this way 1147 1 . To eliminate this random phase, we 
can take the absolute value of the Fourier transform before av- 
eraging. Thus we obtain an effective correlation function: 



C cff ( S ) = (/[0,fc o ]} = < 



(9) 



Note that although it has similar behaviour, our effective cor- 
relation function is expected to be quantitatively different from 
C^ 1 ' (s). Taking the absolute value of the transform reduces the 
cancelling effect between the random thermal phase shifts in 
the statistical average, so that the effective correlation function 
C off (s) decays more slowly with s than C^(s), as calculated 
in the next section. 



4.3 Simulation 

The 1st order correlation function CW (g) can be calculated an- 
alytically 1 33 1, using the theory of |27 1 to account for the phase 
fluctuations. However, it is not possible to obtain an analytic 
expression for the effective correlation function C cS (s) (equa- 
tion l|9}) which we measure. Therefore, to calculate C cS (s), we 
first simulate the quasi-condensate phase fluctuations following 
the theory in |27|. The phase operator is given by: 



<£th (r) = [4ra(r)] 



- 1/2 E 

3 



(10) 



where hj is the annihilation operator of the excitation with 
quantum number j. The solution of the Bogoliubov-de Gennes 
equations for "low energy" excitations (with energies t v < 
hu!±) gives the wavefunctions of these modes: 



/ 7 + (r) = 



'(j' + 2)(2j 



4tt(j + l)R 2 Le 



3)Wr) p (!, 1)(z/i); 



(ID 



where Pj^ are Jacobi polynomials, g = AmK l a/m is the in- 
teraction strength, a is the scattering length, R and L are the 
size of the condensate in the trap and ej = hw z + 3)/4 
is the energy of mode j |40|. We assume the Thomas-Fermi 




Fig. 4. Top: Lines are simulated effective correlation function C c (s) 
(described in the text) as a function of s/L for different T/T,p. From 
the top, T/T<f, = 0,1,3, 10. With increasing T/T^, the curve changes 
from a quasi-gaussian to an exponential and the width decreases. In- 
set: Comparison of C eff (s) (solid line) and C^'(s) (dashed line) at 
T/T^ = 10. Bottom: Width at 1/e of C cB (s) (solid line) and C (1) (s) 
(dashed line) as a function of T/T^. 



approximation for the density n(r), taking R and L from fits 
to images after expansion. To simulate numerically the phase 
fluctuations, we replace the operators a,j and a\ by complex 
Gaussian random variables aj and a* 



,j 3 

These variables have 



a mean value of zero and the correlation (aja*,) = Sjj>Nj 

where Nj = ksT/ (Huj z ( j + 3)/4) is the occupation num- 
ber for the quasiparticle mode j at a given temperature T. We 
assume that the phase fluctuations do not evolve on the time 
scale of the expansion |33|. We have verified this by studying 
images of condensates after the same time-of-flight, but with- 
out the interferometer pulses. In this case, we observe a smooth 
density profile, with no extra features appearing in the Fourier 
transform. 

For a given T/T^, 20 condensates are generated at each 
value of s, with s ranging from 0.2L to 2L. Each condensate 
is integrated over y as in equation l|5}. These simulated absorp- 
tion images are analysed in exactly the same way as the real 
experimental images. The absolute values of the Fourier trans- 
forms of the images are averaged, and the contrast extracted as 
in equation (|9}. The points in Figure [3] show the typical con- 
trast extracted from our simulations for T/T^ — 1. We found 
that these points are very well fitted by the product of a gaus- 
sian and an exponential, for all T/T^. We use these fits as our 
theoretical effective correlation functions C cS (s). 

The effective correlation function was simulated for < 
T/T < 20. Figure presents results for different T/Tj. At 
T = 0, C cff (s) coincides with C^(s) (33). This function is 
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Fig. 5. Example of two experimental effective correlation profiles 
(points) as a function of s/L showing clearly the change in shape 
and width with T/T^, Each point is averaged over 5 condensates, 
each with 2 output port images. Top: At T/T^ = 1.35, the effec- 
tive correlation curve is fitted by a gaussian (solid line), with 1/e 
width Lc/L = 0.62. The fit to an exponential (dashed line) is shown 
for comparison. Bottom: At T /T$ = 4.86, the effective correlation 
curve is fitted better by an exponential (solid line), with 1/e width 
Lc/L = 0.46. The fit to a gaussian (dashed line) is shown for com- 
parison. 

simply the integrated overlap function between the two con- 
densates, and is approximately a gaussian function of the sepa- 
ration s 1321 . As T/Tj, increases, the width of the function de- 
creases and its form gradually becomes exponential. In Figure 
[4j> we plot the 1/e widths of the simulated C cff (s) functions as 
a function of T /T^. For comparison we show also the width of 
(s) 1331 which decreases much faster with T/T^. 



5 Experimental results 

Figure [5] shows two examples of effective correlation curves 
measured using our interferometer and analysed as described 
above. The points shown in Figure |5^ were obtained using a 
magnetic trap with an aspect ratio of 45, and with an atom 
number and temperature corresponding to T/T^ — 1.35, that 
is for small-amplitude phase fluctuations. The points in Fig- 
ure |5j3 were obtained using a trapping aspect ratio of 100, and 
with T/Tcf, = 4.86. The contrast is plotted as function of s/L, 
obtained from a fit to a truncated parabola. Each point corre- 
sponds to an average over 5 condensates. The difference in the 
range of s /L explored is due to different expansion dynamics 
after release from the two different traps. In the more tightly 
confined trap (Figure[5]3), the axial expansion is much slower, 
and thus the fringe spacing decreases more slowly with s. At 
the smallest values of s, it is therefore impossible to measure 
the contrast reliably since we do not observe enough fringes. 
We can extract information about the phase fluctuations from 
both the shape and the width of these effective correlation func- 
tions. 




0.3 b 1 1 1 1 

2 4 6 

Fig. 6. Lc/L as a function of T/T$. The solid curve is a fit using 
the function described in the text. Empty circles correspond to the two 
experimental correlation profiles of Figure|5| Inset: The solid curve is 
the simulated effective correlation function of Figure|4|3. 

5.1 Shape of the effective correlation functions 

First, we observe qualitatively that the shape of the effective 
correlation functions C (s) changes as T/T^ increases. For 
small T/T^, as in Figure^, the curves are clearly gaussian, as 
shown by the fit in the figure. As we increase T/T^, the pro- 
files become rapidly exponential. We see in Figure [SJ) that at 
T/Tfi = 4.86, the curve is already better fitted by an exponen- 
tial than a gaussian. At intermediate values of T/T^, we can 
use the product of a gaussian and an exponential to fit a smooth 
curve through the data. The contribution of the exponential in- 
creases rapidly in importance at finite T, in agreement with 
the simulation, reflecting the increasing amplitude of the phase 
fluctuations with T/T$. 



5.2 Comparison of coherence length L c 

In order to extract quantitative information from the effective 
correlation functions, we define a coherence length Lc, equal 
to the 1/e width of the effective correlation curve C (s). We 
then use this parameter to compare the widths of the measured 
and simulated effective correlation functions. A smooth curve 
is fitted through the data (using the product of a gaussian and 
an exponential) and the 1 /e width extracted. Although the ther- 
mal cloud plays no role in the interference pattern we observe, 
it appears behind the condensates in the p = output of the in- 
terferometer and thus reduces the measured contrast. Indepen- 
dent measurements of the thermal fraction (between 60 % and 
80 % for this experiment) allow us to renormalize the exper- 
imental effective correlation functions to take account of this 
effect. When fitting the renormalized curves, we then fix the 
value at s/L = to unity. 

In Figure[6] we plot Lq/L as a function of T/T^. Impor- 
tantly, we see that the coherence length varies smoothly as a 
function of T/T^, even when the temperature is close to T^, 
This is what we should expect since is simply a charac- 
teristic temperature, defined as that at which the mean square 
fluctuations of the phase difference between two points sepa- 
rated by a distance L is equal to 1 . Therefore it should be borne 
in mind that even condensates at temperatures below are not 
necessarily fully coherent. 

The inset of Figure [6] compares the measured coherence 
widths with the results of the simulation (Figure |4ji). Despite 
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the offset between the two curves, the trend of the data follows 
very well that of the simulation. In Figure |6] the experimental 
data is fitted by a curve A + B exp (—bT/T<fr), giving b = 0.35 
and B = 0.32. This is in reasonable agreement with the simu- 
lation, for which the fit yields b = 0.32 and B = 0.39. In the 
following discussion we consider possible explanations for the 
observed reduction in contrast. 



5.3 Discussion 

As shown in the inset of Figure [6] the measured coherence 
length is offset by about 20% from the results of the simula- 
tion at T/T^ = 0. In order to eliminate various possible causes 
of this discrepancy, and to understand better the limitations of 
our experiment, we have performed several tests. 

Since the accuracy of our experiment relies on the com- 
parison of fringe contrast at different spatial frequencies, it is 
important to take great care in setting up and characterising 
the imaging system. Therefore we measured the modulation 
transfer function (MTF) (see e.g. H48I - ) of our complete imag- 
ing system in situ, using a USAF1940 resolution target en- 
graved with 3-bar square wave patterns of spatial frequency 
4-200 lines/mm, covering the range of spatial frequencies ob- 
served in our interference experiment. By Fourier transform- 
ing images of different regions of the target, we were able to 
compare the magnitude of the different Fourier components in 
the image with those of the target pattern, thereby obtaining 
the MTF of our system. We found that the MTF is approxi- 
mately linear, falling from 1 at zero spatial frequency to at 
118 lines/mm. This resolution limit at 8.5 /im is in agreement 
with earlier characterisations of the system [42 1. The shape of 
the MTF is due almost entirely to the CCD camera |49|, and 
is surprisingly significant at fringe spacings much greater than 
the effective pixel size of 2.5 fim. All contrast measurements 
were corrected by this MTF. In fact, for the data obtained using 
the second trapping aspect ratio, the axial expansion (and thus 
the phase difference developed) was sufficiently small that the 
maximum observed fringe spatial frequency was 16 lines/mm 
and therefore the MTF correction had a negligible effect on the 
experimental effective correlation functions. In the first set of 
data, where the maximum fringe frequency was 38 lines/mm (a 
fringe spacing of 25 /im), the correction was more significant, 
changing the width of the curves by typically 10%, though still 
leaving a 20% discrepancy with the simulation. 

We also considered the error introduced by a small focus- 
sing error. The imaging system is focussed onto the conden- 
sate to within ±0.2 mm by minimising the imaged size of a 
small condensate as a function of the objective lens position. 
However, by varying the time-of-flight used, we may have in- 
troduced an error of up to ±lmm. Such an error could also 
be introduced by small changes to the residual magnetic fields, 
which lead to changes in the release velocity of the condensate 
when the magnetic trap is switched off. We measured effective 
correlation curves for different foci of the imaging system, but 
found that the width of the effective correlation curve changed 
by less than the existing spread in the points. 

Other possible sources of error, such as the alignment of 
the imaging beam with respect to the condensates' fringes and 
correct background subtraction have also been eliminated. 



There remains in our experiment an unexplained pheno- 
menon regarding the distribution of the fringe pattern. We ex- 
pect to see fringes only in the region where the two condensates 
overlap. However, it can be seen in the images of Figure|2]that 
the fringes extend to the edges of each condensate. Moreover, 
these "extra" fringes have the same spatial frequency and phase 
as the central fringes. Although it is possible that a small frac- 
tion of the thermal cloud is coupled by the Bragg beams, inter- 
ference fringes produced in this way 1 30 1 would have a much 
smaller fringe spacing, less than 6 ^m. More importantly, the 
contributions from different parts of the original thermal cloud, 
whose width is ~ 100 [im, would sum incoherently to wash out 
the fringe pattern. It is more likely that these extra fringes arise 
as a result of interactions during the application of the Bragg 
pulses, but better modelling is still needed before evaluating 
whether their presence should increase or decrease the overall 
measured contrast. 



6 Conclusion 

We have demonstrated a new type of matter-wave interferom- 
etry using Bragg beam-splitter pulses and Fourier space anal- 
ysis. Our results show that the expected shape of the correlation 
functions changes from a gaussian-like shape to an exponential- 
like shape when the amplitude of phase fluctuations is increased. 
The coherence length of elongated condensates varies smoothly 
at temperatures close to T&, as predicted by theory. This high- 
lights the fact that the characteristic phase temperature T$ does 
not indicate a transition to full phase coherence, but rather that 
condensates exhibit phase fluctuations at all finite temperatures, 
albeit of small amplitude. This may place constraints on the 
trapping geometries which can be used for creating measure- 
ment devices based on the phase coherence of condensates. 



7 Acknowledgements 

We would like to acknowledge support from IXSEA-OCEANO 
(M.H.), from the BEC2000+ programme and the Marie Curie 
Fellowships programme of the European Union (J.R.), from 
IXCORE (EG.) and from the Chateaubriand Program, CNRS, 
and NSERC (J.H.T.). This work was supported by Delegation 
Generale de 1' Armement, the European Union (grants IST-2001- 
38863 and MRTN-CT-2003-505032) and INTAS (Contract No. 
211-855). 



References 

1. M.H. Anderson et al, Science 269,198 (1995); K.B. Davis et al, 
Phys. Rev. Lett. 75, 3969 (1995); C.A. Sackett and R.G. Hulet, Phys. 
Rev. Lett. 78, 985 (1997); C.C. Bradley et al, Phys. Rev. Lett. 75, 
1687 (1995) 

2. J. Reichel, W. Hansel and T.W. Hansch, Phys. Rev. Lett. 83, 3398 
(1999) 

3. D. Mutter, D.Z. Anderson, R.J. Grow, P.D.D. Schwindt and E.A. 
Cornell, Phys. Rev. Lett. 83, 5194 (1999) 

4. N.H. Dekker et al, Phys. Rev. Lett. 84, 1 124 (2000) 



8 



M. Hugbart et ah; Coherence length of an elongated condensate: a study by matter- wave interferometry 



5. A. Haase, D. Cassettari, B. Hessmo and J. Schmiedmayer, Phys. 
Rev. A 64, 043405 (2001) 

6. J. Fortagh, H. Ott, G. Schlotterbeck and C. Zimmermann, Appl. 
Phys. Lett. 81, 1 146 (2002) 

7. M.P.A. Jones, C.J. Vale, D. Sahagun, B.V. Hall and EA. Hinds, 
Phys. Rev. Lett. 91, 080401 (2003) 

8. W. Hansel, P. Hommelhoff, T.W. Hansch and J. Reichel, Nature 
413, 498 (2001) 

9. H. Ott, J. Fortagh, G. Schlotterbeck, A. Grossmann and C. Zim- 
mermann, Phys. Rev. Lett. 87, 230401 (2001) 

10. T.L. Gustavson et al, Phys. Rev. Lett. 88, 020401 (2002) 
U.S. Schneider et al, Phys. Rev. A 67, 023612 (2003) 

12. Y. Lin, I. Teper, C. Chin and V. Vuletic, Phys. Rev. Lett. 92, 
050404 (2004) 

13. J. Esteve, C. Aussibal, T. Schumm, C. Figl, D. Mailly, I. Bou- 
choule, C.I. Westbrook and A. Aspect, Phys. Rev. A 70, 043629 
(2004) 

14. C. Vale, B. Upcroft, M.J. Davis, N.R. Heckenberg and H. 
Rubinsztein-Dunlop, J. Phys. B 37, 2959 (2004) 

15. EA. Hinds, C.J. Vale and M.G. Boshier, Phys. Rev. Lett. 86, 1462 
(2001) 

16. W. Hansel, J. Reichel, P.Hommelhoff and T.W. Hansch, Phys. 
Rev. A 64, 063607 (2001) 

17. E. Andersson, T. Calarco, R. Folman, M. Andersson, B. Hessmo 
and J. Schmiedmaye r, Phys. Re v. Lett. 88, 100401 (2002) 

18. Y.-J. Wang et al, cond-mat/0407689 

19. S. Dettmer et al., Phys. Rev. Lett. 87, 160406 (2001) 

20. I. Shvarchuck, Ch. Buggle, D.S. Petrov, K. Dieckmann, M. 
Zielonkowski, M. Kemmann, T.G. Tiecke, W. von Klitzing, G.V. 
Shlyapnikov and J.T.M. Walraven, Phys. Rev. Lett 89, 270404 
(2002) 

21. S. Richard, F. Gerbier, J.H. Thywissen, M. Hugbart, P. Bouyer 
and A. Aspect, Phys. Rev. Lett. 91, 010405 (2003) 

22. E.W. Hagley, L. Deng, M. Kozuma, M. Trippenbach, YB. Band, 
M. Edwards, M. Doery, PS. Julienne, K. Helmerson, S.L. Rolston 
and W.D. Phillips, Phys. Rev. Lett. 83, 3112 (1999) 

23. J. Stenger et al, Phys. Rev. Lett. 82, 4569 (1999) 

24. I. Bloch, T.W. Hansch and T. Esslinger, Nature 403, 166 (2000) 

25. D.S. Petrov, G.V. Shlyapnikov and J.T.M. Walraven, Phys. Rev. 
Lett. 85, 3745 (2000) 

26. Y. Castin, R. Dum, E. Mandonnet, A. Minguzzi, I. Carusotto, 
Journal of Modern Optics 47, 2671 (2000) 

27. D.S. Petrov, G.V. Shlyapnikov and J.T.M. Walraven, Phys. Rev. 
Lett. 87, 050404 (2001) 

28. D. Hellweg, L. Cacciapuoti, M. Kottke, T. Shulte, K. Sengstock, 
W. Ertmer and J.J. Arlt, Phys. Rev. Lett. 91, 010406 (2003) 

29. J.E. Simsarian, J. Denschlag, M. Edwards, C.W. Clark, L. Deng, 
E.W. Hagley, K. Helmerson, S.L. Rolston and W.D. Phillips, Phys. 
Rev. Lett. 85, 2040 (2000) 

30. D.E. Miller, J.R. Angl in, J.R. Abo-Shaeer, K. Xu, J.K. Chin and 
W. Ketterle, cond-mat/04 12672 

31. Y. Castin and R. Dum, Phys. Rev. Lett. 77, 5315-5319 (1996) 

32. F. Zambelli, L. Pitaevskii, D. M. Stamper-Kurn and S. Stringari, 
Phys. Rev. A 61, 063608 (2000) 

33. F. Gerbier, J.H. Thywissen, S. Richard, M. Hugbart, P. Bouyer 
and A. Aspect, Phys. Rev. A 67, 051602 (2003) 

34. AA. Michelson, Astroph. J. 51, 257 (1920) 

35. A. Labeyrie, in Progress in Optics XIV, edited by E. Wolf (North- 
Holland, 1976) p. 47-87 

36. Note that the Bragg momentum spectroscopy method does not 
suffer from this problem. 

37. R. Hanbury Brown and R.Q. Twiss, Nature (London) 177, 27 
(1956); 178, 1046 (1956) 



38. B. Desruelle et al, Phys. Rev. A 60, R1759 (1999) 

39. P. Bouyer, A. Aspect, M. Lecrivain, B. Desruelle and V. Boyer, 
patent 00 02704 (2000) 

40. S. Stringari, Phys. Rev. A 58, 2385 (1998) 

41. P. Bouyer, J.H. Thywissen, F. Gerbier, M. Hugbart, S. Richard, J. 
Retter and A. Aspect, J. Phys. IV France 116, 219 (2004) 

42. S. Richard, Ph.D. thesis, University of Paris XI, 2003 

43. F. Gerbier, Ann. Phys. Fr. 29, (2004) 

44. F. Gerbier et al, Phys. Rev. Lett. 92, 030405 (2004) 

45. F. Gerbier, J.H. Thywissen, S. Richard, M. Hugbart, P. Bouyer, 
and A. Aspect, Phys. Rev. A 70, 013607 (2004) 

46. Our cell allows maximum 25 ms free-fall. We increase the time- 
of-flight by applying, after the Bragg pulses, a magnetic field gradi- 
ent opposed to the gravity field. 

47. This would be equivalent to averaging the original images: the 
fringe contrast would be averaged out to the non-interfering density 
profile. 

48. J.W. Goodman Introduction to Fourier Optics 2nd edition, 
McGraw-Hill, 1996. 

49. This has been confirmed by MTF measurements of the Pixelfly 
HiRes camera alone, provided by PCO (www.pco.de). 



